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Convergence of approximate solutions 
of conservation laws 


Sebastian Noelle 1 and Michael Westdickenberg 2 * 


Abstract I 

n this paper we consider convergence of approximate solutions of conser¬ 
vation laws. We start with an overview over the historical developments since 
the 1950s, and the analytical tools used in this context. Then we present some 
of our own results on the convergence of numerical approximations, discuss 
recent related work and open problems. 


1 Exact and approximate solutions of conservation laws 

Hyperbolic conservation laws are used widely to model the local conserva¬ 
tion of physical quantities like density, momentum and energy. Applications 
include gasdynamics, ideal magnetohydrodynamics, shallow water and traffic 
flows. The equations take the form 

du 

— + V • f(u) = 0, for (t, x) G [0, T\ x (1) 

where u(t,x) := ... ,u m (t,x)) is the m-vector of conservative vari¬ 

ables with initial data zt(0, •) =: zto, and f(u) = (/i(u),..., fd(u )) is the flux 
function. It is well known that for nonlinear fluxes, even for smooth initial 
data, the solution of (1) may cease to exist in the classical sense due to the 
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formation of shocks in finite time (cf. Riemann 1859 [43]). Therefore it is 
necessary to consider weak solutions, i.e. functions u that are bounded and 
satisfy (1) in the distributional sense. Existence of weak solutions for general 
systems of conservation laws is largely unknown, and weak solutions are not 
unique. However, for scalar equations (with m = 1) the Cauchy problem is 
well understood. 

As usual, a possible strategy to establish existence of weak solutions of the 
conservation law (1) is first to regularize the problem and obtain a sequence 
of approximate solutions {f( £ } e for e > 0, then to show that a subsequence 
of u e converges to a limit function u which is a weak solution of (1). The 
regularization mostly used is the vanishing viscosity method, where u e solves 
the parabolic problem d t u e + V • f{u e ) = eAu e . Since in general exact solu¬ 
tions of conservation laws are not known, approximations through numerical 
computations are of utmost importance in the applications. As the grid size 
tends to zero, again a sequence of approximate solutions is created, and then 
a convergence analysis is important to make sure that the results are reliable. 
We say more on that issue in Section 3. 

We will assume that the sequence {u e } e of approximate solutions is uni¬ 
formly bounded in L°°([0, T] x R d ). For the vanishing viscosity method and 
for many numerical schemes this can be obtained from a (discrete) maximum 
principle or an invariant region argument. Then it is possible to extract a sub¬ 
sequence which converges weak* to some limit function u € L°°([0, T] x R d ). 
But in order to prove that u is a weak solution of (1), it is necessary to show 
that the weak limit of f(u e ) coincides with f(u). For a weakly converging se¬ 
quence u e and a nonlinear flux / this is not true in general because oscillations 
may occur. A sufficient condition, however, would be strong L^-convergence. 
So the question is how this additional information can be obtained. 

We already mentioned that for a general nonlinear hyperbolic conserva¬ 
tion law there may exist many different weak solutions corresponding to the 
same initial data. In order to ensure uniqueness one imposes an additional 
condition, often called an entropy condition because of its analogy with the 
second law of thermodynamics, which selects the physically relevant weak 
solution out of all possible ones. For scalar conservation laws several related 
formulations are in use. We give the one that relies on convex entropies. 

We call a pair of functions rj: R m —> R and q: R m —> R d an entropy-entropy 
flux pair if r] is convex and if the compatibility relations 

q' k {u) = rj\u)f' k {u) for all u G R m and k = 1,..., d (2) 

hold. Then the entropy condition reads as follows: Out of all weak solutions 
of the conservation law (1) given above, we only select those satisfying 


dr](u) 

dt 


+ V • q(u) < 0 


in distributional sense 


( 3 ) 
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for all entropy-entropy flux pairs ( r],q ). We also assume that the initial en¬ 
tropy is bounded: f Rd r/(uo) dx < oo. These solutions are called weak entropy 
solutions. For scalar conservation laws all convex functions can serve as 
entropies, and it turns out that weak entropy solutions are unique. This is 
exceptional. For systems of conservation laws there often exists only a lim¬ 
ited number or even only one single (the physical) entropy, because the cor¬ 
responding compatibility relations are much more restrictive. The entropy 
inequalities provide crucial estimates in the convergence analyis for {u e } e , 
and the lack of enough entropy-entropy flux pairs is the main reason why the 
existence theory for systems of conservation laws is not as well developed as 
for the scalar case. For more on the entropy condition we refer to Lax [30] and 
Liu [34]. Further information on hyperbolic conservation laws can be found 
in Lax [29], Smoller [46] and Dafermos [12]. 

Our paper is organized as follows. In Section 2, we give an overview over 
the developments in the theory of conservation laws since the 1950s, and 
the analytical tools used to study the convergence of approximate solutions. 
These include classical compactness arguments, approximation theory, weak 
convergence methods (compensated compactness and measure-valued solu¬ 
tions), and the kinetic approach. In Section 3, we present some of our own 
work on the convergence of numerical approximations. Finally, in Section 4, 
we discuss recent related work and open problems. 


2 Historical remarks on compactness arguments 

In this section we give an overview over some analytical tools used to study 
convergence of sequences of approximate solutions for conservation laws. 


2.1 Regularity estimates and compactness 

The classical approach to prove strong L i 1 oc -convergence of a sequence {u e } e 
of functions is to prove that the u e possess some positive regularity, i.e. the 
sequence is uniformly bounded in some appropriate function space. If the 
function space is BV 1 the space of functions with bounded variation, then 
strong compactness in Lj follows from Helly’s Theorem. 

In 1957 Ladyzhenskaya [27] and Oleinik [39], both former students of 
Petrovskii in Moscow, independently published results on a scalar conser¬ 
vation law with strictly convex flux in one space dimension. The entropy 
condition they used to select the physical solution (today known as “condi¬ 
tion E ”) can be interpreted as a one-sided Lipschitz condition: For any time 
t > 0, the solution u(t, •) is Lipschitz continuous where it is increasing (in x). 
There is no restriction where it is decreasing. In particular, decreasing jump 
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discontinuities are allowed. From this a BV bound can be derived. We would 
like to remark that the one-sided Lipscliitz continuity was already observed 
by E. Hopf [20] in his pioneering paper on the Burgers equation, but he did 
not use it systematically as selection criterion. While Ladyzhenskaya and 
Oleinik both studied the vanishing viscosity limit, Oleinik together with her 
student Vvedenskaya [54, 39] also proved convergence of the Lax-Friedrichs 
scheme [28]. In 1963, Oleinik’s paper was translated into English, while La¬ 
dyzhenskaya and Vvedenskaya’s work was not. Oleinik’s contribution is also 
discussed in full detail in Snroller’s famous book [46]. This may explain why 
condition E is nowadays usually attributed to Oleinik. 

The multidimensional scalar problem was studied by Conway and Smollcr 
[9] and Vol’pert [53] in the 1960s, and by Kruzkov [25] in 1970 who proved 
existence and uniqueness of weak entropy solutions for the scalar conservation 
law. Kruzkov’s main observation was that the solution operator of the scalar 
conservation law is an Id-contraction: He used the special convex entropies 

r)(u, k ) := | u — fc| and q(u, k ) := sgn(w — k)\f(u) — f(k)\ (4) 

(nowadays known as Kruzkov entropies), where fc £ I is a real parameter, 
and a “doubling of variables”-argument to conclude that 

/ \u(t,x) — v(t,x)\dx < / |tto(a;) — Vo(x)\dx (5) 

Jm d Js. d 

for all times t > 0. Here u and v are weak entropy solutions of (1) corre¬ 
sponding to initial data uo and vq, respectively. Uniqueness of weak entropy 
solutions of (1) is an immediate consequence of this estimate. And since the 
solution operator of (1) is invariant under translations in x, we can apply (5) 
also to v(t, x) := u{t,x + h) for any h £ K. d and obtain from this a HV-bound 
|| u(t, -)l|w(R d ) < ll u o||w(R‘ i )- An analogous estimate for temporal regularity 
can then be derived using (1). So Helly’s theorem can be applied. 

Let us briefly mention that there is also a well-developed theory for sys¬ 
tems of conservation laws in one space dimension relying on BV bounds. This 
theory started with the work of Glimnr [17] in 1965 who introduced a nu¬ 
merical scheme that produces approximate solutions with bounded variation. 
Later on, Bressan and his coworkers extended Glimm’s ideas considerably and 
set up a very powerful framework for studying existence and uniqueness of 
systems of conservation laws in 1-d. We refer to Bressan [3]. 

BV bounds also play a major role in the design of higher order accurate 
numerical algorithms. In 1983, Harten [19] introduced the class of TVD (to¬ 
tal variation diminishing) schemes and proved their convergence in the linear 
case. Oslier [40] introduced the class of “E-schemes”, the most general class of 
schemes which satisfy a discrete entropy inequality, and proved convergence 
of second order accurate semidiscrete TVD schemes for nonlinear scalar con¬ 
servation laws. Convergence of a second order accurate fully discrete scheme 
was proved by Oslier and Tadmor in 1988 [41]. 
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2.2 Approximation theory 


Approximation theory in the context of scalar conservation laws means to 
estimate the distance between an approximate solution (at some time t > 0, 
say) and the uniquely defined weak entropy solution which is already known 
to exist. This approach gives error estimates and convergence rates. Applied 
to numerical schemes it tells you how close an approximate solution is to the 
exact one. 

The first who gave error estimates for approximate solutions of scalar con¬ 
servation laws in several space dimensions was Kuznetsov [26], who studied 
the vanishing viscosity method and the Lax-Friedrichs scheme. His analysis 
relies on the ZA-contraction property of the solution operator of (1) we dis¬ 
cussed in the previous section. In fact, it is possible to obtain an estimate 
like (5) even if v is not an exact, but only a (suitable) approximate solution. 
Then additional terms measuring the approximation error enter on the right 
hand side which can be controlled if the initial data is in BV. More precisely, 
one can estimate the ZA-distance between the weak entropy solution u and 
the approximate solution u e at a given time t in terms of the distance and 
HV-norm of the initial data, the residual etc. These estimates depend on 
e (which could be the gridsize in numerical schemes, for example), and the 
convergence rate proved by Kuznetsov is e 1 / 2 . 

Kuznetsov’s ideas were used later on to give convergence rate estimates for 
more sophisticated numerical schemes, as well. We refer to Sanders [45], Vila 
[52], Cockburn, Coquel and LeFlocli [7], Noelle [37], and Cockburn, Gremaud 
and Yang [8]. We also refer to Bouchut and Perthame [1] who reformulated 
Kruzkov and Kuznetsov’s approximation theory in a general, versatile form. 

We also mention that recently, Tadinor [50] and coworkers developed a 
somewhat different approach to proving error estimates and convergence rates 
for scalar conservation laws in 1-d with strictly convex flux, the so-called 
Liphtheory. The idea is to measure the distance between the exact and the 
approximate solution not in the ZA-norm, but in the much weaker topology 
of Lip^K.). This is the topological dual of the space of Lipschitz continuous 
functions. The approach can also give pointwise error estimates: The theory 
shows that the convergence rate at some given point x £ R only depends on 
the regularity of the exact solution u in a small neighborhood around x. 


2.3 Weak convergence methods 

For modern numerical schemes using unstructured grids it may be too hard 
or even impossible to prove positive regularity in terms of boundedness of a 
sequence {u e } e in appropriate function spaces. So the classical approach of 
proving compactness in Lj oc by using Helly’s theorem or Sobolev embeddings 
breaks down. There are other tools, which run under the name of Weak 
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Convergence Methods, that can prove strong compactness without showing 
regularity first. These methods are specialized in the sense that they rely 
heavily on the structure of the problem at hand. 

In 1954, Lax [28] proved compactness in L 1 of the solution operator of a 
one-dimensional, not necessarily strictly convex, scalar conservation law us¬ 
ing a weak topology. More than twenty years later, Tartar [51] and Murat 
[35] introduced an even more powerful method, applicable also for systems of 
conservation laws, the compensated compactness theory. A typical result is 
the div-curl-Lemma: One assumes that two sequences of vector-valued func¬ 
tions {U e } e and {V e } e are bounded in (Lj ! oc (R ti )) jv and that the divergence 
of U e and the curl of V e are both precompact in Hr(M. n ). Then there exist 
subsequences such that 

U e ■ V e —> U ■ V in distributional sense. ( 6 ) 

Note that the assumptions on the differentiability of U e and V e are not 
sufficient to obtain strong compactness in L^ oc (R") (from which ( 6 ) would 
follow trivially). But in combination they give enough information to have 
weak continuity of the scalar product. 

The div-curl-Lennna has been used by Tartar to study the convergence of 
sequences of approximate solutions {w e } e of the scalar conservation law ( 1 ) 
in one space dimension. Compensated Compactness arguments can also be 
used to study certain systems of conservation laws in 1-d. As an example, we 
refer to the work of DiPerna [13] and Chen [5] on isentropic Euler equations. 

Another weak convergence method is DiPerna’s theory of measure-valued 
solutions for multidimensional scalar conservation laws. As mentioned before, 
we assume that the sequence {n e }e of approximate solutions is uniformly 
bounded in L°° and extract a subsequence converging weak* in L°°. Then 
the sequence converges strongly in Lj oc if and only if no oscillations occur, i.e. 
if there is convergence pointwise a.e. The concept of measure-valued solution 
was first introduced by Tartar [51] in 1975, who used Young measures to 
describe the oscillations occuring in {n e } e . A Young measure associated to 
{u e } e is a weakly measurable mapping v from [0,T] x K. d into the space of 
probability measures Prob(R), such that for all continuous functions / we 
have f{u e ) —Kv,f> in distributional sense, where the pairing <v,f> is 
given by 

<v,f>{t,x) := [ f(\)dv( ttX) {\). 

J R 

So the distributional limit of any nonlinear function of u e can be written down 
using one single u. A Young measure is called a measure-valued solution of 
the Cauchy problem (1) if the following identity holds 

d 

— < v, id > + V• < v, / >= 0 in distributional sense. (7) 

Consistency with the entropy condition then means: For all convex 77 
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— <v,p> +V- <v,q>< 0 in distributional sense, (8) 

where q is the corresponding entropy flux. 

Note that a measure-valued solution is even weaker than a weak (distribu¬ 
tional) solution. One can show that for any bounded sequence of approximate 
solutions {u e } e there exists a Young measure v and a subsequence converging 
to it in an appropriate sense. That v is a measure-valued solution of (1) then 
follows from the consistency of the approximation. 

One has <v, id>= u, where u is the distributional limit of (it £ } e . But in 
general <v,f>^ /(u). Equality holds for arbitrary nonlinear / if and only 
if the Young measure reduces to a Dirac measure, i.e. if v is of the form 
u (t,x) — 3u(t,x) with S u the Dirac distribution on R centered at u. This is 
equivalent to saying that the sequence {w £ } € converges strongly in Lj oc . 

The following result is due to DiPerna [14] (see also Szepessy [47]) 

Theorem 1. Assume that initial data uq € L 1 0 L 00 (R <i ) is given and that 
there is exists a Young measure v satisfying the following properties 

(i) the function (t,x) lX ), |id|> is in L°°([0, T\, ld(R d )), 

(ii) v is a measure-valued solution of (1), 

(in) v is consistent with the entropy condition, 

(iv) v assumes the initial data uq in the following sense 

lim - f f <v^ s x ), |id — uo^)! > dsdx = 0. (9) 

t_i-0 t J Q J R d 

Then the Young measure reduces to a Dirac measure, i.e. V(t, x ) = &u(t,x ) f or 
a.e. (t,x) € [0,T] x R d , where u is the unique entropy solution of (1). 

Assumption (iv) says that the initial data must be attained in a stronger 
sense than just the sense of distributions. It excludes the occurence of oscilla¬ 
tions in the sequence of approximate initial data {uo} e - Then condition (iii) 
assures that no oscillations can develop at later times: the only oscillations 
that can exist in the sequence (u e } e are those transported into the system 
from the initial data. If the initial data converge strongly, so does {u € } e at 
any later time. 

Theorem 1 was used to prove convergence of the streamline diffusion shock¬ 
capturing method by Szepessy [48], finite difference methods by Coquel and 
LeFloclr [10], spectral viscosity approximations by Chen, Du and Tadmor [6] 
and finite volume schemes on unstructured polygonal grids by Kroner and 
Rokyta, [24], Kroner, Noelle and Rokyta [22], Noclle [36], and others. Rohde 
[44] extended the method to weakly coupled systems of conservation laws. 
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2.4 Kinetic formulation and velocity averaging 


The kinetic formulation was introduced by Lions, Perthanre, and Tadmor 
[32]. They show that there is a one-to-one correspondence between weak 
entropy solutions of the scalar conservation law (1) and solutions of a linear 
transport equation with source term, for which a certain nonlinear constraint 
holds. More precisely, one considers a “density-like” function p depending on 
(t,x) € [0, i] x R d , and on an additional variable net, which is a solution 
of the transport equation 


dp 

dt 


+ f'(v) • Vp 


dm 

dv 


in distributional sense. 


( 10 ) 


Here rn is a nonnegative bounded measure. Equation (10) is supplemented 
with an assumption on the structure of p. If y: R —> R is defined by 


x( v \ a ) 


+1 if 0 < v < a 
— 1 if a < < 0 

0 otherwise 


for a £ R, then p should have the form 


P(t, x,v) = x(v\u(t,x)) 


( 11 ) 


( 12 ) 


for some scalar function u. If p is a, solution of the kinetic equation (10) and 
satisfies the nonlinear constraint (12), then the function u is the unique weak 
entropy solution of (1). Vice versa, if u is the unique entropy solution of (1), 
then there exists a nonnegative bounded measure m such that the function 
p defined by (12) is a solution of the kinetic equation (10). 

One big advantage of the kinetic formulation for scalar conservation laws 
is the possibility to apply velocity averaging lemmas to obtain regularity and 
compactness for (sequences of approximate) weak entropy solutions of (1). 
One can show that the macroscopic quantity u has more regularity than the p 
whose ^-average it is, see Golse, Lions, Perthanre and Sentis [18] and DiPerna, 
Lions and Meyer [15]. 

In Westdickenberg and Noelle [55] we used the kinetic formulation together 
with the velocity averaging technique to prove the convergence of a class of 
finite volume schemes for scalar conservation laws in several space dimensions. 
We will discuss this result in more detail in Section 3.3 below. 

We remark that already in 1984, Brenier [2] constructed and analyzed an 
approximate evolution operator for (1) based on a kinetic decomposition of u. 
In 1998 Perthanre [42] used the kinetic formulation to study again uniqueness 
of weak entropy solutions and error estimates in the sense of Kuznetsov’s 
approximation theory. It is also possible to give a kinetic formulation for 
certain systems of conservation laws, e.g. the isentropic Euler equations in 
1-d. In Lions, Perthanre, Tadnror [33] and Lions, Perthanre, Souganidis [31] 
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this kinetic formulation was used together with compensated compactness 
to prove existence of weak entropy solutions, thereby extending DiPerna [13] 
and Chen’s [5] results mentioned above. 


3 Some convergence results for finite volume schemes 


In this section we give an overview over some of our own results on the conver¬ 
gence of finite volume schemes. These schemes define approximate solutions 
which are piecewise polynomial on a given polygonal grid of and in gen¬ 
eral discontinuous at the cell interfaces. Integrating the conservation law (1) 
over a cell K with faces e C dK from time t n to time t n+1 = t n + At gives 
the following update formula for the cell average u : 


u 


n+1 

K 


= U 


K 


At 

W\ 


E M 9 U K,e- 

eC dK 


(13) 


Here g 7 ^ e is a numerical flux which is consistent with /(it) • riK,e , where 
riK, e is the outward unit normal of face e. This flux has to be conservative, 
satisfy some upwinding property (e.g. an .E-flux), and it may be higher or¬ 
der accurate in space and time, see e.g. [22, 36]. We consider explicit time 
discretizations, which means that the fluxes g^ can be computed directly 
from the cell averages u'^, in a neighborhood of K. Such discretizations are 
only stable for timesteps At which are so small that waves originating from 
face e = K D K', where the approximation is discontinuous, do not cross the 
neighboring cells K and K' completely during one timestep. If one imposes 
this CFL-condition (named after Courant, Friedrichs and Lewy [11]), then it 
is possible to prove a discrete maximum principle. 


3.1 Discrete entropy inequalities 


The discrete maximum principle is enough to assure that a subsequence con¬ 
verges weak* to some limit function, but it does not guarantee that this limit 
is a weak solution of the conservation law. The crucial additional estimate 
needed to prove convergence is a discrete entropy inequality of the form 


V(u 


r>-•>(<*) + W\ 


£ 

eO dK 


| G\ e < Ch 2a . 


(14) 


Here h is the maximal diameter of a cell of the computational grid, a € (^, 1] 
is a constant depending on the triangulation and the polynomial reconstruc¬ 
tion, the constant C depends on the entropy rj and its derivatives, and G\ e 
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is a numerical entropy flux which is consistent with the entropy flux q ■ riK,e 
and is closely related to the numerical flux g^ e . 

In [22] we proved an entropy inequality of this type for a class of higher 
order schemes based on the Lax-Friedrichs and the Engquist-Osher schemes, 
and in [36] we could treat schemes which extend Godunov’s exact Riemann 
solver to higher order accuracy (see [22, 36] for the definitions and the original 
references). 

There are several points in the derivation of estimate (14) which are worth 
recalling. In [36] we first rewrite the update (13) as a convex combination 
of one-dimensional contributions. This is natural, since the numerical fluxes 

e approximate one-dimensional fluxes / • riK,e propagating in the normal 
direction to the cell-faces. Adapting the weights of the decomposition to the 
local wave speeds, our analysis admits larger time steps than previous work, 
even for the classical first order schemes in 1-d [49]. 

Next we use Tadmor’s decomposition [49] of any E-flux as a convex com¬ 
bination of Lax-Friedrichs and Godunov’s fluxes. If one uses a piecewise con¬ 
stant reconstruction, corresponding to a first order accurate discretization, 
then the discrete entropy inequality with (7 = 0 follows naturally from the 
upwinding properties and the CFL-condition. For higher order polynomial 
reconstructions, we have to generalize Harten’s TVD conditions discussed in 
Section 2.1 above to multidimensional reconstructions. The estimate for the 
higher order Lax-Friedrichs numerical flux follows readily, but for the higher 
order Godunov’s flux, which yields the best resolution of shock discontinu¬ 
ities, the structure of the Riemann solution has to be exploited in full detail, 
see [36]. Once the one-dimensional discrete entropy inequalities are obtained, 
convexity arguments yield the multidimensional inequality (14). 


3.2 Convergence and error estimates via 
L 1 -contraction 


Equipped with the discrete maximum principle and entropy inequality (14), 
we can now prove convergence and error estimates. In [22] we gave the first 
prove of convergence for higher order accurate finite volume schemes on un¬ 
structured grids for initial data which are merely in L°°, using DiPerna’s 
theory of measure-valued solutions. Subsequently, Kroner and collaborators 
designed practical algorithms along the lines of our convergence results. 

In closely related work, Cockburn, Coquel and LeFloch [7] used Kuznet¬ 
sov’s approximation theory to prove convergence and error estimates for BV 
initial data in 1994. In [37], we could generalize these error estimates, and 
simplify their proof, using our discrete entropy inequalities. 

Even though inequality (14) is the crucial estimate, careful work is required 
to finish the convergence proof. In this process, many papers (including [22]) 
assume that the triangulation is regular, meaning that the ratio of the outer 
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diameter to the inner diameter of a cell is bounded above. In [36], we removed 
that restriction and treated grids where the cells may become flat in the 
limit as h —> 0. In [37], we quantified the influence of irregular grids on the 
convergence rate. An example in [36] shows that the rate with which the grid 
may degenerate is optimal. 


3.3 Convergence via kinetic formulation 

In Westdickenberg and Noelle [55] we gave a new convergence proof for se¬ 
quences of approximate solutions of the scalar conservation law, produced 
by higher order finite volume schemes on unstructured grids in several space 
dimensions. We did not use the theory of measure-valued solutions, but relied 
on the kinetic formulation and velocity averaging lemmas instead. The key 
tools in our analysis were once more the discrete entropy condition (14) and 
the following compactness theorem. 

Theorem 2. Let 1 < p < 2 and 0 < 7 < 1. Fix a compact set K Ct and 
assume that there are sequences {p e } e , {?n e } e; and {7T e } e uniformly bounded 
in ZA(R d+1 x and L 1 (K,C*{R d+1 )), resp., such that 

d p € dm e 

+ f'( v ) ' V/j £ = —- h 7 r £ in distributional sense (15) 

((t,x) £ M d+1 , v £ K). If the following nondegeneracy condition holds 

sup nreasju £ A: r + f'(v) ■ f = 0} = 0, (16) 

(r,{)eR‘ i + 1 

then the sequence {z e } e belongs to a compact subset of Lj oc (R d+1 ), where 

z e ■= f p e (-,v)dv. (17) 

JK 

Here, Al(R d+1 ) is the space of bounded Radon measures, and C'*(R d+1 ) is the 
topological dual of the space of Holder continuous functions. This theorem is 
a variant of the compactness theorem in [32]. Assumption (16) means that 
the flux function must be nonlinear enough. The regularizing effect stated in 
Theorem 2 does not exist for linearly degenerate problems (e.g. for advection 
equations). Note that the m e on the right hand side of (15) is a bounded 
measure, as we should expect from the kinetic formulation (10). The quantity 
7r e is an error term (e.g. the numerical residual). It is measured in a function 
space with negative regularity, hence in a rather weak topology. 

To prove Theorem 2 we decompose z e into two parts and show that one 
part can be made arbitrarily small in Lj using the nondegeneracy of the flux 
(16), while the other part has some positive regularity in terms of Sobolev 
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(in fact Besov) spaces and is therefore strongly compact in L] oc . Applying 
Theorem 2 to numerical approximations we obtain the strong compactness 
needed to pass to the limits in nonlinear quantities as explained in Section 1. 


4 Related work and open problems 

In [23, 38] Kroner and Ohlberger derived a posteriori error estimates based 
on Kuznetsov’s approximation theory and recent related work of Gallouet, 
Herbin, Clrainais-Hillairet et al. [4]. Ohlberger [38] developed a fully adap¬ 
tive, implicit finite volume scheme for scalar convection-reaction-diffusion 
equations based on these estimates. 

In a recent preprint, Hwang and Tzavaras [21] used the kinetic formulation 
together with the velocity averaging technique to study the convergence of 
approximate solutions of scalar conservation laws that are produced by either 
a relaxation method or by a diffusion-dispersion approximation. 

To conclude this note, we would like to sketch some important open prob¬ 
lems concerning the convergence of numerical approximations. The conver¬ 
gence rate of ft, 1 / 2 proved by Kuznetsov for the Lax-Friedriclis scheme is only 
optimal for first order approximations of discontinuous solutions of linear 
advection equations. The rate obtained so far for unstructured grids is only 
ft 1 / 4 . A first attempt to overcome this barrier may be found in [8]. For strictly 
convex scalar conservation laws one expects a rate of ft log ft, but this conjec¬ 
ture is so far only backed up by numerical experiments. For schemes which 
are formally higher order accurate one expects higher convergence rates away 
from discontinuities. For strictly convex scalar conservation laws a result in 
this direction is proven in [16]. It is also shown there that for systems of con¬ 
servation laws, characteristics of one family crossing a numerical shock layer 
belonging to a different family may carry first order errors into the smooth 
postshock region, so the convergence rate deteriorates there even for higher 
order schemes. In ongoing work G. Kreiss et al. are beginning to design special 
numerical schemes which attempt to filter out this pollution effect. 
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